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In this brief technical communication it is demonstrated how using Bethe ansatz technique the 
explicit expression for the two-point free energy distribution function in (1+1) directed polymers 
can be derived in rather simple way. Obtained result is equivalent to the one derived earlier by 
Prolhac and Spohn 0]. 
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INTRODUCTION 



We consider the model of directed polymers described in terms of an elastic string cj)(r) directed along the r-axes 
within an interval [0,£] and defined by the Hamiltonian 



g; H[<Kt),V\ = J rfr{i[a r 0(r)] 2 + l/[^(r),T]}; (1) 



where the disorder potential V\<j>, r] is Gaussian distributed with a zero mean V(<j>,t) = and 



V((j>, t)V{4>', t') = u5{t - t')5{4> - 4>') (2) 

The parameter u describes the strength of the disorder. The partition function with the fixed boundary condition, 

4>{i) = x is: 



Z(,-,= / * Vl'ie '"^ = i'S V i-.m.r)) 13) 

o 



ZM = / T)Mt\ P-PHM = 

'0(O)=O 



Correspondingly, F(x) is the free energy of the polymer which at time t arrives to the point x. In the limit of large t 
£N| ' random free energy scales as 

>'• 

^H ■ pF = (3f t + Xf (4) 

r-- ; 

where /q is the trivial self-averaging contribution (which can be easily eliminated by simple redefinition of the total 
^■^ ' free energy and the partition function) and 

O; A = ^« J t) 1/8 «* 1/8 (5) 

In the limit t — > cc the random quantity / ~ 1 in eq.(Q| is described by the universal Tracy- Widom distribution 04a] ■ 
. — • The aim of the present brief communication is the study the two-point free energy probability distribution function: 

W(f 1 J 2 ;x 1 ,x 2 ) = lim Prob[/(zi) > A; /(a*) > h] (6) 

Some time ago the result for this function has been derived in terms of the Bethe ansatz replica technique under a 
particular decoupling assumption [lj. Here I'm going to recompute this function, again in terms of the same general 
scheme of the Bethe ansatz approach but using somewhat different computational tricks (which do not require any 
supplementary assumptions). Since this function depends only on the distance between the two points, x = \x2 — #i|, 
to simplify formulas I'll consider the particular case: x\ = —\x and x 2 = +hx- In other words, instead of ([5]), I'll 
concentrate on the probability distribution function defined as follows: 

W(h,fcx) - lim Prob[/(-x/2) > / i; /(as/2) > f 2 ] (7) 

Just recently it has been proven Q that the result of the present calculations (see eqs. (|52")) -(l5!J ]) below) is equivalent 
to the ones obtained earlier [l(. 



II. TWO POINT DISTRIBUTION FUNCTION 

In terms of the partition function, eq.®, above the probability distribution function, eq.([7]) can be defined as 
follows: 



(—l) L (— l) H f itr -,R 

Wtfufcx) = lim V V l l cxp(A£/ 1 + XRf 2 ) Z(-x/2) exp{(3f t} Z(x/2) exp{/?/ i} (8) 

A— >oo L — ' L — ' Li\ it: L J L 

L=0 R=0 

Here, the averaging, denoted by (...) is performed over random potentials ([2]). Performing the standard averaging of 
the (L + i?)-th power of the partition function, eq.©, one gets 



W{h,h;x)= lim J2J2 



(-1) L (-1) 



R 



L=0 R=Q 



- exp(Ai/i + \Rf 2 + /3(L + R)f t) ^(-x/2, ..., -x/2, x/2, ..., x/2; t) (9) 



L R 

where the time dependent wave function ^(xi, ..., xn] t) is the solution of the imaginary time Schrodinger equation 

N , N 



W(x;t) = [^ £#. +5 «$>(*« -a%)] *(*!*) 



o=l 



a^b 



where k = /3 3 u and the initial condition 



tf(x;t = 0) = II? =1 ($(a; a 



(10) 



(11) 



A generic eigenstate of such system is characterized by N momenta {q a } (a = 1, ...,JV) which split into M (1 < 
M < N) "clusters" described by continuous real momenta q a (a = 1,...,M) and having n a discrete imaginary 
" components" 



q a = q" = Qa- y( n " + 1 - 2r) 



(r = l,...,n a ) 



with the global constraint 



M 



J2 n « = N 



(12) 



(13) 



A generic time dependent solution ^(x, t) of the Schrodinger equation (|10[) with the initial conditions, eq. (jll[) . can 
be represented in the form of the linear combination of the eigenfunctions ^q (x): 



N i r r 
% %;«)= J2w p (M) (q,* 



M—\ 



where we have introduced the notation 



|C M (q,n)| 2 vI/W(x)* q M )*(0) exp{-£ M (q,n)*} 



/ M f r +oo t » 1 M 
„ 1 «/ — OO „ i i 



(14) 



(15) 



;( M ), 



and 8(k, m) is the Kronecker symbol. For a given set of integers {M; n\, ...., tim} the eigenfunctions "Jq (x) can be 
represented as follows: 



JV 



^ m) (x)=eii 



V a<b 



1 + in 



sgn(z a - x b ) 
1v a - qr b 



exp 






(16) 



where the summation goes over JV! permutations V of JV momenta q a , eg. (1121) . over JV particles x a ; the normalization 
factor 



,N 



|C M (q,n)| 2 



M 

n 



f(n a -ri/j)! 



W! ria = l( Kn «) a</3 1 0a 



; K + 71/3)1 



(17) 



and the eigenvalues: 



E 



M 



M i 2 



24,3 



-iV 



(18) 



The last term in the above expression is the self-averaging part of the free energy; choosing /o = k 2 /(24/3 2 ) this term 

drops out of the further calculations. Note also that according to the definition, cq. (|16[) . ^q (0) = N\ 
Substituting eas.([I3 1) -(fTg j) into eq.©, we get: 



W(hJ 2 ;x) 



lim > 



(-1) 



L+R 



A^oo ^ L\R\ 

L+R>1 



Q\Lh+\Rf 2 



(19) 



L+_R 



1 



M 



M=l Q = l 



E 



+ 00 



_dqa_ 
2Trnn a 



; n a g-|n Q ?» + 54jti' 



A/ 



<5(Vn Q , i + i?) |C M (q,n)| 2 



a=l 



>< EEEnn 

■p{L,R) -p(L) -p(R) a =l c =l 



<7„(i) — gw-R) — »« 



fj,W — 9-p(H) 



exp 



a— 1 c— 1 



where 



I Cm (q,n) 



-r-r \q a - qp - —{n a - np)\ 
a <p \la -qp - l -j{n a +np)\ 



(20) 



In ea. (|19p the summation over all permutations V of (L + R) momenta {qi, ..., Ql+r} over L "left" particles {xi, ..., xl} 
and R "right" particles {yR, ■•■, yi} split into three parts: the permutations V^ L > of L momenta (taken at random out 
of the total list {gi, ...^l+r}) over L "left" particles, the permutations J>( R ) of the remaining R momenta over R 
"right" particles, and finally the permutations V^ L ' R ^ (or the exchange) of the momenta between the group "L" and 
the group " R" . It is evident that due to the symmetry of the expression in eq. (|19|) with respect to the permutations 
■p( L ) and V^ R ' the summations over these permutations give just the factor L\ R\. 

Further simplification comes from the following general property of the Bethe ansatz wave function, eq. (fl"6|) . It has 
such structure that for ordered particles positions (e.g. x\ < x 2 < ... < xjv) in the summation over permutations the 
momenta q a belonging to the same cluster also remain ordered. In other words, if we consider the momenta, eq.(|12p. of 
a cluster a, {qf , gf , ..., g" o }, belonging to the particles {xi t < Xi 2 < ... < Xi n }, the permutation of any two momenta 
q" and q® of this ordered set gives zero contribution. Thus, in order to perform the summation over the permutations 
*p(.L,R) - m eQ pT)| [i [ s sufficient to split the momenta of each cluster into two parts: {gf,...,g^ \\q^ +1 .,.,q™ }, where 
m a = 0, 1, ..., n a and where the momenta q" , ..., q^ belong to the particles of the sector "L" (whose coordinates are 
all equal to —x/2), while the momenta q^ +1 ...,q% belong to the particles of the sector "R" (whose coordinates are 
all equal to +x/2). 

Let us introduce the numbering of the momenta of the sector " R" in the reversed order: 



->• 



-> 



1 l 
9 2 



where m c 



and (si. eq.flU])) 



1m a + l 



->• 



(21) 



q*" = q a + -z-(n a + 1 - 2r) = q a + —{m a + s a + l- 2r) 

By definition, the integer parameters {m a } and {s a } fulfill the global constrains 

M 



^2m a = L 

=i 

M 

y^gq = R 



M 



(22) 

(23) 
(24) 



In this way the summation over permutations J>( L ' R ^ in ea. (TT9t is changed by the summations over the integer 
parameters {m a } and {s Q }, which allows to lift the summations over L, R, and {n a }. Straightforward calculations 
result in the following expression: 



M 



00 (—iW 



M\ 

M=l Q = l 



00 - H 

_ ^1 j— 



dq a 



.m Q +s Q >1 



( ^ r ^ r * i 1 * / \ 2 K * / S3"! 

expj Am a /i + As Q / 2 - -xm a q a + -xs a q a - -«m a s„ - ■77n\ m a + s Q )g a + ^T7(m« + s Q j j 



27tk(to q + s a ) 

24/T 



x |CM(q, m + s)| 2 G M (q, m, s) 



(25) 



where 



GM(q,m, s) = nnn v.;. *nnn 



a— 1 r— 1 r' — 1 
M 



9? " q** 



a</3r=l r' = l 



q'r - <l*r> ~ iK 
. 8? " <l*r> , 



T-r T(l + m a + s a ) 

jj r(i + m Q )r(i + SQ ) x 



>< n 

a</3 



r 


-, , rn a +mf}+s a +sp , i /•„ 

1 t 2 "r k vy Q 


9/3) 


r 


, — rria+mp+Sa 

1 ' 2 


^ + *(fe- 


-9/3) 


r 


1 T 2 "+" k V^" 


-9/3 


y 


r 


1 ' 2 


^ + ^(fe" 


-10) 



After rescaling 



<fa -> Yx 1 a 



2A 2 
a; — > x 

K 



1/^V3 = l ( ^ 2<)1 /3 



with 



2 V 
the normalization factor |C*M(q, m + s)| 2 , eq. (l2T))) . can be represented as follows: 

,~ 2 T-r I A (m a + s a ) - A (mfl + as) - «p Q + «P/3 1 
|C M (qm + s)| = [[L- - - ■ — U = 

a</3 |A(to q + s a ) + \{mp + sp) - ip a +ipp\ 



= f ff [2A(m Q + s a )] ) x det — 



t) - Wa + \(mp + Sp) + ipp 
Substituting eqs. (j2"Tl) - (|30l) into eq. (j2"5")) and using the Airy function relation 



= 1,...,M 



exp 



; A 3 (m Q + s Q )^ 



dy Ai(j/) exp X(m a + s a )y 



we get 



W( fl , h ; X )= Urn 1 + £ i-J- fl 



A— >oo 



M! 

M=l a=l 



+00 



dy a dq a 



2vr 



Ai(y Q + (£) E (- 1 



m a +s a — 1 



m Q +s a >l 



expJAm Q (y„ + /1 - -Jfeai) + Xs a (y a + f 2 + -ifea;) - \ 2 m a s a x\ 



(26) 



(27) 
(28) 

(29) 



(30) 



(31) 



x det if 



(\m a , Xs a , q a ); (Xmp, Xsp, q fj ) 



a,0=l,...,M 



Gm 12A' 



(32) 



where 



K 



(Am, As, g); (Am , As , q ) 



(33) 



Am + As — iq + Am' + As' + iq 1 
Using the relation 

exp{-A 2 msa;} = / ' g /2 3 expj--^ - -gf - ^3 + Am\/z£i + As\/^6 + iK m + a)v^6J (34) 

the expression in eq.([32|) can be represented as follows: 



00 (—-\) M 

w{h,h-,x)= 1 +e 4?r n 

M=l ' a=l 



dy Q dg Q d$,i a d& a d£ 3a 



2ir (2tt) 3 / 2 

x exp|--(£i Q + 2^\/^) 2 ~ 2 fe Q ~ 2 iqa ^ 2 ~ 2^ 3 "J 
where 
5Af(q,y,li,li,/i,/ 2 ,a;) = lim J| 



SM(q.,y,$i,ti,h,f2,x) (35) 



A^oo 



x detX 



J^ (_!)»»-+«- 1 exp{Am Q ( 2 / Q + / 1 + ei a ^) + As Q (y Q + / 2 +6 Q \/^)} 

,m Q +s Q >1 

I (Am Q , As Q , g Q ); (Am^j, As,g, ^) 



a,/3=l,...,M 



'M 



/«q 

l2A' m ' 



(36) 



To demonstrate how the summations over {m Q } and {s a } are performed in the limit A — > 00 let us consider the 
example of a general type: 



M 



R = lim TT 



E (-1)"" 1 exp{An a y a } 



$(A; ni, ..., n M ) 



(37) 



where $ is a function which depends both on A and on all summation parameters {n\ , ..., Um}- The above summations 
can be represented in terms of the integrals in the complex plane: 



M 



R = lim F[ 



1 f dz a 



2i J c sin(7rz Q ) 



exp{\z a y a } 



$(A; zi, ..., z M ) 



(38) 



where the integration goes over the contour C shown in Fig. 1(a). Shifting the contour to the position C shown in 
Fig. 1(b) (assuming that there is no contribution from infinity), and redefining z a — > z a /X, in the limit A — > 00 we get: 



A I 



R=U 



I &Z(x r i 

/ exp{z Q y Q } 

JC z a 



1 f dz, 

2m j L 



lim $(A; zi/A, ..., z M /X) 

A— ^-00 



(39) 



where the parameters y a and z a remain finite in the limit A — > oo. 

The double summations over m a and s a in eq. (|36|) can be represented as follows 



00 00 



00 00 



E (-i)™-- 4 -- 1 = E *-.o E (-ir*- 1 + E <w ° E (-i) 5 ^ 1 - E (- 1 )" 1 "" 1 E (- 1 )'" -1 ( 4 °) 

'm Q +s Q >l s«=0 m a — l m a —0 s Q — 1 m a — l s a — 1 

Thus in the integral representation, eqs. (j3"7|) - (f39|) , for the function in eq . ([36]) . we get 

M 



<5jvf(q,y,$ 1 ,$i,/i,/2,x) = J]. 



Q = l 



dz\ a dz2 a (2m' 



(2m\, . 27T? ., . 1 \ 

I o(z 2q ) H o(z ia ) I 

\z la z 2a z la z 2a J 



c ( 27Ti ) 2 ^ z i. 
exp|zi Q (2/ Q + /i +£i a V%) + z 2a {y a + h +& a Vx)\ 



x 



x dct K[(z la , z 2a ,q a ); (z 1/3 , Z2j3,q{s)] x 



lim Gm 

A— >og 



/Kq zi 
V2A' A' 



(41) 



(42) 



c 



12 3 4 5 



12 3 4 5 



(a) 



(b) 



FIG. 1: The contours of integration in the complex plane used for summing the series: (a) the original contour C; (b) the 
deformed contour C'; 



Using the explicit form of the factor Gj./, eq. (|26p . and taking into account the gamma function property limi z i_>o T(l- 
z) = 1, we find 






(43) 



Thus, in the limit A — > oo the expression for the probability distribution function, eq.(|35]), takes the form of the 
Fredholm determinant 



W(f x ,fr,x) = 1 +J2 



(-1) 



A I 



M 

Ml H 

M=\ Q = l 



x dy a dq a d£ la d& a d& a , 2 ,-, 



2tt 



(2tt) 3 /2 



exp { _ 2^ lQ + 2* 9a ^ 2 ~ 2^ 2q _ 2^ Q ^ 2 ~ 2^ 3 "} 



dz la dz 2a (2m-. 2m 1 
, I <5(z 2 J + S(z la ) I 

V ( 2m ) Vz 1q Z2q «l «2ct' 



x exp|zi Q (y Q +/i +6aV^) +^2 Q (2/ Q +/2 + 6 a V^)| 

x det K[(z la , z 2a , q a ); (z 1/3 ,z 2/3 , qp)] a/3=1 M 

= det[l - A] 



with the kernel 



A[(z 1: z 2 , q);(zi, z 2 ', <?')] = 



+oo 



dy d£id£2d£ 3 



Ai(y + q 2 - i&Vx) 



2tt (2tt) 3 / 2 

cxp {~2 (^ + 2 iq ^ 2 ~ 2^ 2 ~ 2 iq ^ 2 ~ 2^ 3 / 

1 



'2m' 



( i-Kl 27TI 

< I S(z 2 ) + 6{zi 



2m 

z 2 



z\z 2 



x exp 



{zi (y + h + &VS) + z 2 (y + f 2 + £ 2A /i) } 



1 



Zi + z 2 - iq + Zi + z 2 + iq 
In the exponential representation of this determinant we get 



W(h,f 2 ,x) = exp[- ]T i 

M=l 



TrA 



a/ 



(44) 



(45) 



(46) 



where 



TrA 



M 



M 

n 



Ai(y + q 2 a -i&Vx) x 



+ °° dydg a d£id&d£3 
-oo 2n (2tt)3/2 

1/. 1. «2 1- 1 . ^ N 2 I 



exp|--(6 + -iq a Vx) 2 - 2(6 - ^^V^) 2 ~ 2^ 3 J 

dzi Q dz 2Q /2?ri 2?ri 1 \ 

d(z 2a ) H o(«ia) x 

\Zl™ Zt„ Zi„Zo„/ 



la ( 27 ™) 2 ^ z la 
x exp|zi a (y + /i +£iy/x) + z 2a (y + f 2 + &Vx)j 



1 



M 

J-J- zi a + z 2a -iq a + zi a+1 + z 2a+1 +iq a +i 
Here, by definition, it is assumed that Zi M+1 = Zi t (i = 1,2) and <?a/+i = q\- Substituting 
1 



z \a + z 2q — *<Za + z l«+l + ^2q + 1 + ^Qa+1 

into eci.(j4"T|). we obtain 



doj a exp 



[Zla + z 2 a ~ iq a + Z la+1 + Z 2a+ i + iq a +l) U a 



r°° M 

TtA M = / dW\...dw M TT U.(w a ,Wa+l 



(47) 



(48) 



(49) 



where 



A(uj,uj') 



+00 



dydq d£ 1 d£ 2 d£3 



Ai(y + g 2 + cj + u' - i&y/x) x 



2tt (2tt) 3 / 2 

exp \ 2^ + 2 iq ^ 2 ~ 2^ 2 ~ 2* 9 ^ 2 ~ 2^ 2 ~ ^ W ~ ^J 

dz\dz 2 (2-ki 2iri 1 

d(z 2 ) H <5(zi) 

V zi Z9 Zi 2 



^122 



l c , (2ni) 2 V zi 
x expjzi (y + fi+ £iVx) + z 2 (y + f 2 + ^v 7 ^) } 
Integrating over z\ and z 2 we get: 

A(w,w / ) 



(50) 



+ OO Z' + OO 



^ / o~~ Ai (y + "? 2 ~~ /1 + w + w ' + 9 ^) exp i ~^( u - UJ ')} 





+00 /> + oo 

dy 

-OO /"+00 



2tt 



2~ Ai(y + <? 2 - / 2 + w + J - -iqx) exp{ -ig(w - u/) } + 



' 2tt 



OO ./ —OO 



+ °° d^d&d^ , 2 , 

/ 9 W2 Ai(j/ + q +w + w -i&V^J x 
-00 l^ 71 "/ 



exp|--(6 + -jiqVxf - -(& - -jiqVx) 2 ~ 2"6 2 - «?(w - w')J 



0(y + A + 6Vi) 0(j/ + / 2 + 6 >/i) 



(51) 



where 0(y) is the step function. Redefining, £1 = (t — rj)/>/2, £ 2 = (t + r])/y/2, £3 = (it + £)/\/2, and integrating 



over g, t and £, we find the following result: 

A(u,,«/) = 2 1 /3x[2 1 /3( w _/ 1 ) ; 2 1 /3( w '_/ 1 )] cxp {I( CJ _ CJ ') a; } + 



2 1 '*K 



2 1 / 8 ( W -/a) s 2 1 /8( w '-/ 2 )]e J q>{-i(a;-a;')a J } 



r+oo 

2 2 / 3 / dy 



+ OG 



c/?7 



2tt 



Ai 



'"(, + „-^)]« i ^( f+ -' + ,J5) 



{-I" 2 - r» - ?^ W >4©>M- V!M 9+/2+ VI) 



x exp 



where /i, 2 = |(/i,2 ~ x 2 /16) and i^(w, w') = Jp 00 dy Ai(y + ui) Ai(y + u/) is the Airy kernel. 
Thus the distribution function W(/i, f2)x), eq.([7|), is given the Fredholm determinant 

W(h,f 2 -x) = det[l - i] 

where A is the integral operator with the kernel A(u,uj') (uj,uj' > 0) given in eq.([52|). 

Note that using explicit expression ([52]) one can easily test the obtained result for three limit cases: 



lim A{uj,{jj') 



2 l ' 3 K 



V3( w _ / 2 ) , 2 1 /3( w ' _ / 2 )j exp{-i( W - «')*} 
V3( w _ £) , 2V8( W ' - A)] exp{i( w - «>} 
Km%.') - 2^ 2 1 />-/ 1 /2), 2 1 / 3 (^'-/ 1 /2)l 0(/i - / 2 ) + 



lim A(u,w') = 2 1/3 K 



x-yO 



2 1 /3^r 2 V3( w _ /2/2 ) , 2 V3( W ' _ / 2 /2)l 0(/ 2 - A) 



(52) 



(53) 



(54) 
(55) 

(56) 



which demonstrate that in the case /i i2 — > — oo we recover the usual GUE Tracy- Widom distribution for / 2i i cor- 
respondingly, while in the limit case x — > we find the usual GUE Tracy- Widom distribution for f± (in the case 
h > fz) and for / 2 (in the case / 2 > /i), as it should be. 

III. CONCLUSIONS 

In view of the recent proof [6j that the result of the present calculations is equivalent to the one obtained earlier 
by Prolhac and Spohn |l| we can conclude that the Bethe ansatz replica technique has demonstrated (once again) 
its the efficiency and robustness which allows to perform computations of sufficiently complicated objects in rather 
simple way. 
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